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The Chinese rufous horseshoe bat ( Rhinolophus sinicus) has been suggested to carry the direct 
ancestor of severe acute respiratory syndrome (SARS) coronavirus (SCoV), and the diversity of 
SARS-like CoVs (SLCoV) within this Rhinolophus species is therefore worth investigating. Here, 
we demonstrate the remarkable diversity of SLCoVs in R. sinicus and identify a strain with the 
same pattern of phylogenetic incongruence (i.e. an indication of recombination) as reported 
previously in another SLCoV strain. Moreover, this strain possesses a distinctive 579 nt deletion 
in the nsp3 region that was also found in a human SCoV from the late-phase epidemic. 
Phylogenetic analysis of the Orfl region suggested that the human SCoVs are phylogenetically 
closer to SLCoVs in R. sinicus than to SLCoVs in other Rhinolophus species. These findings 
Received 3 September 2009 reveal a closer evolutionary linkage between SCoV in humans and SLCoVs in R. sinicus, defining 

Accepted 11 December 2009 the scope of surveillance to search for the direct ancestor of human SCoVs. 


Severe acute respiratory syndrome (SARS) is an infectious 
disease in humans caused by a coronavirus (CoV) named 
SARS CoV (SCoV) (Rota et al, 2003). Small carnivores, 
including civet cats and raccoon dogs, were found to carry 
SCoVs and were therefore believed to be the immediate 
sources of the SARS epidemic in 2003 (Guan et al., 2003; 
Kan et al, 2005). Later, Rhinolophus species (horseshoe 
bats) were found to carry a diverse group of SARS CoV-like 
CoVs (SLCoV) and, thus, may be the natural reservoir of 
SLCoVs and SCoVs (Lau et al., 2005; Li et al., 2005). 
However, the currently sampled SLCoV in bats (Bt- 
SLCoV) may not be the direct ancestor of SCoVs in 
humans (Hu-SCoV) (and thus the civet and raccoon dog), 
because Bt-SLCoVs are phylogenetically distant from Hu- 
SCoVs (Ren et al, 2006). Nonetheless, our previous study 


The GenBank/EMBL/DDBJ accession numbers for the sequences 
obtained in this study are FJ588686 (complete genome of Rs672), 
FJ588687-FJ588692 (nsp3, nsp13, S, Orf3a and N sequences of 
Rs806) and FJ588681-FJ588685 and FJ561 748 {cytB sequences of 
the hosts of Rp3, Rml, Rf 1, Rs672, Rs806 and HKU3). 

tThese authors contributed equally to this paper. 


demonstrated the potential recombinant origin of a Bt- 
SLCoV strain (Rp3) (Hon et al, 2008), and we speculated 
that an uncharacterized SLCoV lineage may exist in its host 
that may represent the direct ancestor of Hu-SCoV. It 
should be noted that we previously misidentified the host 
of Rp3 as Rhinolophus pearsonii (Li et al, 2005); this is now 
corrected to Rhinolophus sinicus based on phylogenetic 
analysis of its cytB sequence in this study (Fig. la). 
Therefore, we anticipate that an investigation of the genetic 
diversity of SLCoVs in R. sinicus could provide insights 
into the possible direct ancestor of Hu-SCoVs. 

Faecal samples of 24 R. sinicus individuals were collected 
from distinct locations within China during September 
2006. RNA was extracted and used for first-strand cDNA 
synthesis as described previously (Li et al, 2005). Viral 
detection using PCR amplification of a 440 nt conserved 
region of nspl2 (RNA-dependent RNA polymerase) was 
performed as described previously (Lau et al, 2005). For 
host species identification of positive samples, the mito¬ 
chondrial cytochrome B gene {cytB) was sequenced as 
described previously (Cui et al, 2007). The cytB genes of 
the hosts of the previously reported Bt-SLCoVs were also 
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Fig. 1 . Phylogeny of the hosts of Bt-SLCoVs and other Rhinolophus species based on their cytB sequences, (a) Bayesian 
phylogeny constructed using mrbayes (Ronquist & Huelsenbeck, 2003) under the best-fit substitution model determined by 
modeltest 3.7 (Posada & Crandall, 1998). The taxa are represented by GenBank accession numbers; the hosts of Bt-SLCoVs 
are highlighted with diamonds. BtCoV279 and BtCoV273 (GenBank accession numbers DQ648857 and DQ648856) were 
excluded, since they were isolated from the same field samples as Rml and Rfl, respectively (S. Zhang, personal 
communication). Bayesian posterior probability support (BPS) of selected nodes is indicated as percentages. Species names 
are indicated to the right. Bar, 0.03 nucleotide substitutions per site. The geographical origins of R. sinicus taxa are indicated by 
circled numbers, referring to the legend in (b). (b) Geographical origins (at the level of province) of the R. sinicus taxa. Black 
lines on the map represent provincial borders. 


sequenced. Two positive samples, designated Rs672 and 
Rs806, were collected from two relatively distant locations, 
in Guizhou and Hubei provinces, respectively (Fig. lb). We 
obtained their genome sequences by using PCR amplifica¬ 
tion and sequencing strategies described previously (primer 
sequences available upon request) (Li et al, 2005). For 
Rs672, the sequences for all putative coding regions were 
obtained. For Rs806, we failed to recover its full-length 
genome sequence, and only the regions covering nsp3, 
nspl3, S, Orf3a and N were sequenced (Fig. 2b). The cytB 
sequences of the hosts of Rs672, Rs806, Rp3, Rml, Rfl (Li 
et al., 2005; Ren et al., 2006) and HKU3 (Lau et al., 2005) 
were analysed with those of other available Rhinolophus 
species. The phylogeny suggests that the hosts were 
correctly identified in the original publications with the 
exception of the host of Rp3 (Li et al, 2005; Ren et al, 
2006), as mentioned above. Notably, the R. sinicus 
sequences appear as two topologically distinct clusters in 
the phylogeny, i.e. the hosts of Rp3, HKU3 and Rs806 in 
clade 1 and the host of Rs672, with several other taxa, in 
clade 2 (Fig. la). The geographical origins of clades 1 and 2 
cover five provinces in southern China (Fig. lb). 

A plot of similarity (Lole et al, 1999) across the genome of 
Rs672 suggests that its 5' region (i.e. Orfla and Orflb ) is 
more similar to that of Hu-SCoV than to that of Bt- 
SLCoV. In contrast, the 3' region of the genome of Rs672 is 
generally more similar to Bt-SLCoV than to Hu-SCoV (Fig. 
2a). This pattern of incongruence is similar to that 


observed in the previously reported potential recombinant 
Rp3 (Hon et al, 2008). We therefore applied the 
recombination detection methods implemented in gard 
(Kosakovsky Pond et al, 2006) and lard (Holmes et al, 
1999) to locate the potential recombination breakpoint in 
Rs672. The parameters and strategies employed are the 
same as described previously (Hon et al, 2008), except that 
Rs672 replaced Rp3 as the query. Both analyses suggested 
that the phylogenetic incongruence in the genome of Rs672 
is statistically significant, and the potential breakpoint was 
estimated at the nucleotide immediately after the start 
codon of the spike gene (model average support in gard 
>0.8; P value in likelihood ratio test in lard <0.0001). 
This potential recombination breakpoint in the genome of 
Rs672 is at the same position as the one detected in the 
genome of Rp3 (Hon et al, 2008). The genome regions 
upstream and downstream of this potential breakpoint 
were respectively designated the major and minor parental 
regions. 

Two phylogenies were constructed based on the coding 
sequences of the parental regions of selected Hu-SCoV 
(«=9) and all available Bt-SLCoV («=6) (Fig. 2c). To 
ensure the quality of the alignment of divergent regions, we 
used concatenated codon alignment instead of full genome 
nucleotide alignment (Fig. 2b). For Rs806, both phyloge¬ 
nies suggest its close phylogenetic relationship with another 
Bt-SLCoV in R. sinicus, HKU3. For Rs672, its phylogenetic 
positions in the two phylogenies are incongruent. Based on 
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Fig. 2. Phylogenetic incongruence within the 
genome of Rs672. (a) Similarity plot based on 
concatenated codon alignment as indicated in 
(b). ‘Hu-SCoVs’ refers to all Hu-SCoVs in (c); 
‘Bt-SLCoVs’ refers to HKU3, Rml and Rf 1. 
The plot was generated using Simplot (Lole 
et al., 1999), with window size of 2000 nt and 
a step size of 200 nt, under the F84 distance 
model, (b) Schematic diagram of the concat¬ 
enated coding sequences used in (a) and (c). 
ORFs are displayed as filled arrows and the 
non-structural protein products of Orfla and 
Orflb are displayed as shaded boxes. For 
annotation of ORFs and non-structural 
proteins, refer to the sequence of TOR2 
(GenBank accession no. AY2741 19). 
Sequenced regions of Rs806 are also indi¬ 
cated. (c) Bayesian phylogenies constructed 
using concatenated codon alignments of the 
indicated regions. Unavailable regions of the 
sequence of Rs806 were treated as missing 
data in the Bayesian estimations. The phylo¬ 
genies were constructed as described in the 
legend of Fig. 1. Hu-SCoV and Bt-SLCoV 
lineages are indicated. Accession numbers, 
host species [civets (C), humans (H) or bats 
(B)] and strain names are given. Bar, 0.01 (left) 
and 0.02 (right) nucleotide substitutions per 
site. Asterisks represent nodes with BPS 
>90%. Phylogenetically incongruent taxa in 
both phylogenies are highlighted by diamonds. 


Fig. 2(c) (left panel), the major parental regions of Rs672 
and Rp3 (designated the HB-SLCoV lineage) shared a 
monophyly that is phylogenetically closer to Hu-SCoVs 
than to Bt-SLCoVs. In contrast, the minor parental region 
of Rs672 is clustered with the Bt-SLCoVs lineage (Fig. 2c, 
right panel). The current data suggest that Rs672 and Rp3 
are likely to have evolved from a common ancestor with 
the same pattern of phylogenetic incongruence. Interest¬ 
ingly, the host of Rs672 is phylogenetically distinct from 
those of the other three Bt-SLCoVs in R. sinicus (Fig. 1), 
implying that Rs672 and Rp3 (as well as their hosts) may 
have diverged a relatively long time ago, which may also 
explain the observation that the minor parental region of 
Rs672 did not share a monophyly with that of Rp3. 

Despite the statistical significance, such incongruence 
could be explained as laboratory artefacts or asymmetrical 
genome evolution and not necessarily as homologous 
recombination. Nonetheless, the identical pattern of 
phylogenetic incongruence identified in Rp3 and Rs672 
may exclude the possibility that it represents a laboratory 
artefact. Alternatively, if the mode of evolution (e.g. 
substitution rate) of Orfl in Rp3/Rs672 is significantly 
different from that of other Bt-SLCoVs, it could also lead 
to an apparently closer phylogenetic relationship between 
the Orfl of Rp3/Rs672 and HuSCoVs. However, from a 


biological perspective, it may not be sensible to assume that 
the same ORFs in viruses of the same host species would 
evolve in a significantly different manner, as Rs806, HKU3, 
Rp3 and Rs672 all reside in the same host species, R. 
sinicus. In other words, unless differential selection 
pressure was operating on the viruses residing in specific 
individuals of R. sinicus, the homologous genome regions 
of Rs806, HKU3, Rp3 and Rs672 would be expected to 
evolve in a similar manner. Therefore, based on the current 
data, homologous recombination seems to be the most 
plausible explanation for the observed phylogenetic incon¬ 
gruence. If such recombination did happen, it is unlikely to 
have occurred recently, since Rp3 and Rs672 showed 
considerable divergence, and their major parent is expected 
to be a Bt-SLCoV that shares a close phylogenetic relation¬ 
ship with Hu-SCoVs. 

We previously used a Bayesian relaxed molecular clock to 
estimate the time of divergence between Rp3 and Hu- 
SCoVs (tDiv-Hu/HB), and the results suggested that the 
interspecies transfer of SLCoVs to the amplifying host (e.g. 
civets) might have happened a median of 4.08 years before 
the SARS outbreak (Hon et al, 2008). Here, we have 
reanalysed the dataset, using the same method and model 
parameters as described previously (Hon et al, 2008), with 
the addition of two sequences derived from this study 
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(Rs672 and Rs807), to investigate whether the addition of 
these highly relevant sequences would alter our previous 
estimate significantly. Fig. 3 shows a time-scaled phylogeny 
constructed using beast version 1.4.8 (Drummond & 
Rambaut, 2007). The time of divergence between Rs672/ 
Rp3 and Hu-SCoVs (i.e. tDiv-Hu/HB) was estimated to be 
at a median of 1999.04, with 95 % highest posterior density 
(HPD) from 2002.06 to 1992.87, which is consistent with 
out previous estimate (median 1998.51, 95% HPD 
1993.55-2001.32) (Hon et al, 2008). The length of branch 
A in Fig. 3, representing the window period between the 
hypothetical interspecies transmission event and the onset 
of the SARS epidemic, was estimated at a median of 4.29 
years (95% HPD 0.70-9.57 years). The estimate in this 
study is fairly consistent with our previous estimate 
(median 4.08 years, 95% 1.45-8.84 years) (Hon et al, 
2008), although its HPD range is slightly wider. This 
analysis was repeated with exponential and logistic 
coalescent tree prior assumptions and the results were 
generally consistent (not shown), suggesting that our 
estimation is robust under different coalescent tree prior 
assumptions. The robustness of the estimation of this 
relatively short window period further supports our 
previous speculation that this uncharacterized HB-SLCoV 
lineage may contain the direct ancestor of SCoV. 

The mp3 region of Orfla is relatively variable in the 
genome of Bt-SLCoVs, as a number of large deletions were 
observed within this region (Ren et al, 2006). Interestingly, 


in PCR amplifications of the mp3 region of Rs672, we 
identified a size variant in one of the amplicons, which was 
later confirmed (using sequencing) to contain a deletion of 
579 nt (Del-579nt). In an attempt to characterize the 
potential viral subpopulations in sample Rs672, we tried to 
culture and isolate viruses using methods described 
previously (Li et al, 2005), but were unsuccessful. 
Nonetheless, our sequencing results suggest that the 
flanking sequences of Del-579nt in both size variants of 
the amplicon are identical (GenBank accession no. 
GQ870287), and these results were confirmed using 
another set of primers for PCR amplification and 
sequencing (not shown). Moreover, the sequences of 
cloned amplicons (n=10) in the S region of Rs672 were 
almost identical, except for some non-recurrent substitu¬ 
tions which are probably due to PCR artefacts or the 
presence of a homogeneous viral quasispecies population. 
Therefore, there is no evidence to support the presence of 
multiple distinct viral subpopulations in sample Rs672, 
except for Del-579nt in the mp3 region. We therefore 
believe that the consensus genome of Rs672 is a robust 
representation of the homogeneous viral population in the 
sample. Surprisingly, Del-579nt was also found in an Hu- 
SCoV strain in the late phase of the 2003 epidemic, 
ShanghaiQXC (GenBank accession no. AY463060), but not 
in other Hu-SCoVs or Bt-SLCoVs available in GenBank. 
Considering the observation that the viral genomes in 
Rs672 are homogeneous except for Del-579nt in the mp3 
region, this deletion was likely to have occurred recently 
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Fig. 3. Estimation of the time of divergence between the major parental region of Rs672/Rp3 and Hu-SCoVs. The time-scaled 
phylogeny was constructed as described previously (Hon et at., 2008), using beast version 1.4.8 (Drummond & Rambaut, 
2007) under the assumption of a constant coalescent tree prior and an exponentially relaxed molecular clock (uced) (Drummond 
et al., 2002, 2006). The dataset was the same as described previously (Hon et al., 2008) except for the addition of the full 
sequence of Rs672 and the partial sequence of Rs806 (indicated by filled circles); the unsequenced region of Rs806 was 
treated as missing data in the Bayesian analysis. tDiv-Hu/HB is indicated by a filled diamond. The dashed line indicates branch 
A, the window between the hypothetical interspecies transmission event and the onset of the SARS epidemic. Shaded bars 
represent 95% HPD of the estimated divergence time of the corresponding nodes (with Bayesian posterior probability >70 %). 
The taxa are annotated as in Fig. 2, except that the sampling time is given at the end. 
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within the host of Rs672; in other words, Del-579nt in 
Rs672 and ShanghaiQXC might have arisen independently 
and was probably not acquired through homologous 
recombination. We cannot establish a plausible hypothesis 
for the origin of Del-579nt; a wider survey of Del-579nt in 
Bt-SLCoVs may shed light on its origin. 

In conclusion, this is the first study to report the 
intraspecies diversity of Bt-SLCoVs. Our data suggest the 
presence of at least two distinct genotypes of Bt-SLCoV in 
R. sinicus (i.e. Rp3/Rs672 and HKU3/Rs806); the diversity 
of Bt-SLCoV in Rhinolophus macrotis and Rhinolophus 
ferrumequinum is currently unknown. Our findings also 
provide evidence for the potential recombinant origin of 
Rp3 and Rs672. Their major parent has a relatively closer 
phylogenetic relationship with Hu-SCoVs, suggesting the 
possible presence of a Bt-SCoV lineage in R. sinicus that 
may contain the direct ancestor of Hu-SCoVs, as we 
proposed previously (Hon et al, 2008). However, we note 
that these speculations are based on a few strains only, and 
extensive surveys on the prevalence of this genotype would 
indicate their credibility. R. sinicus is an extremely 
common Rhinolophus species in China (Lau et al., 2005), 
and our data imply that Hu-SCoVs are phylogenetically 
closer to some Bt-SLCoVs in R. sinicus than to those in R. 
macrotis and R. ferrumequinum. Therefore, we suggest 
more-focused surveillance of Bt-SLCoVs in R. sinicus, 
which may provide insights into the possible direct 
ancestor of Hu-SCoVs. 
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